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We study S-matrix correlations for random matrix ensembles with a Hamiltonian H = Ho + f, 
in which Ho is a deterministic N x N matrix and ip belongs to a Gaussian random matrix ensemble. 
Using Efetov's supersymmetry formalism, we show that in the limit N —* oo correlation functions 
of S-matrix elements are universal on the scale of the local mean level spacing: the dependence of 
^v^j ■ Ho enters into these correlation functions only through the average S-matrix and the average level 

density. This statement applies to each of the three symmetry classes (unitary, orthogonal, and 
symplectic) . 
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I. INTRODUCTION 



The energy levels and/or the scattering matrices of a variety of physical systems with randomness (e.g., complex 
nuclei, disordered conductors, classically chaotic systems, etc.) exhibit universal behavior: the statistical properties 
of the observables can be separated into the universal parts and the non-universal parts specific to individual systems. 
C| There have been increasing evidences that the universal parts depend only on the fundamental symmetries of the 
underlying Hamiltonian and are well described by a random matrix ensemble with a Gaussian distribution (see Ref. (]J 
for a review). According to the fundamental symmetries, there are three classical random matrix ensembles: Systems 
with broken time-reversal symmetry are described by the unitary ensemble and time-reversal invariant systems by 
either the symplectic or the orthogonal ensemble depending on whether spin-orbit coupling is present or not 0. In 
spite of their many successful applications, random matrix models lack firm foundation. Especially, the Gaussian form 
of the probability distribution is used for mathematical convenience and is not motivated by physical principles. It is 
therefore necessary and important to investigate whether statistical properties are identical for more general forms of 
the probability distribution consistent with fundamental symmetries. 

There has been several work along this direction. Hackenbroich and Weidenmuller || considered a non-Gaussian and 
unitary invariant probability distribution: P{H) cx exp [-NtTV(H)], where N is the dimension of the Hamiltonian 
matrix H and V{H) is independent of N and arbitrary provided it confines the spectrum to some finite interval and 
generates a smooth mean level density, in the limit N — > oo. For each of the three symmetry classes, using Efetov's 
\ supersymmetry formalism, they showed that both energy level correlation functions and correlation functions of S- 
matrix elements are independent of P(H) and hence universal if the arguments of the correlators are scaled correctly. 

For realistic situations, it is likely that the Hamiltonian is not completely random but contains some regular parts. 
If the total Hamiltonian H — Hq + p where Hq is a deterministic part and ip is a random one, the probability 
G ' distribution takes a unitary non-invariant form: 

P{H) oc ex$[-NttV{<p)] = exp [-JVtr V{H - H )\ . (1.1) 

For the unitary ensemble, Brezin et al. discussed the universality of 2-point energy level correlations for V(ip) — 
O ■ p 2 /2 + gp 4 . General n-point energy level correlation functions were shown to be universal by Brezin and Hikami || 
for V{p) oc ip 2 . (The other type of unitary non-invariant distrubution P(H) oc exp{— A^tr \V(H) — HHo]} was also 
■ considered by Zinn- Justin ||.) 

' R ecen tly, we M numerically found the same universality of the 5*-matrix correlations for the distribution function 
Eq. (1.1) with Vyp) oc ip 2 for the orthogonal ensemble, i.e., with the average S- matrix S been taken as the parameters, 
the correlations are independent of Hq while S depends on Hq. Our purpose of the present article is to analytically 
show this universality in any of the three symmetry classes. More precisely, we show 



1111 /> <- * [S* cjdj * =Mu,p(E),S(E)), (1.2) 



where m, n, fcj, lj are non negative integers, the bar denotes ensemble average. The universal functions fp depend 
on the symmetry classes (/3 = 1, 2, and 4 for orthogonal, unitary, and symplectic classes) and are independent of Hq, 
except for the indices {a^, hi, hi, Cj, dj, lj}, while the average local level density p and the average ^-matrix S depend 
on Hq. 
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II. THE MODEL 



Following the approach of Ref. ||, we write the scattering matrix S(E) as 



S ab {E) = S ab - 2m^W^ aii \d{E)~ x 



(1 . V 



(2.1a) 



in which 



D{E) = E + i0 + - H + mWW\ 



(2.1b) 



E is the energy, + is positive infinitesimal, H represents the projection of the full Hamiltonian onto the interaction 
region, and W describes the coupling between the eigenstates of the interaction region and the scattering states in 
the free-propagation region. The indices a, b refer to the physical scattering channels, and /i, v refer to the complete 
orthonormal states characterizing the interaction region. 
We assume that N x N matrix H can be written as 



H = H Q + ip, 



(2.2) 



where Hq is a given, nonrandom, Hermitian matrix, and ip is a member of the Gaussian ensemble. The symmetry 
property of Hq is the same as ip. The independent elements of the matrix <p are uncorrelated random variables with 
a Gaussian probability distribution centered at zero. The second moments for the unitary ensemble are given by 



A2 
N 



(2.3) 



(See Ref. M for the orthogonal and the symplectic cases.) Here, A is a strength parameter. 



III. DERIVATION 



For definiteness, we show the derivation for the unitary ensemble and m,n < 2 in Eq. (1.2). The generalization 
to the other symmetry classes and/or higher values of m and n is straightforward and commented upon in Sec. IV. 
The derivation is based on the use of Efetov's supersymmetry method |||l0| . We take the notation from Ref. || and 
use the [1, 2] block notation for the matrix representation in which 1 and 2 refer to the retarded and advanced block, 
respectively. 

Consider the following generating function: 



Z{J) 



det [D p (E p ) + 2nWJ p {F)W^] 
det [Dp(E p ) - 2irWJ p (B)W < <] 



(3.1) 



where D p (E p ) = diag [D{E 1 ),D^{E 2 )] , J P (F) = diag [Ji(F), J 2 (F)], and J p (B) = diag [J X (B), J 2 (B)}. The scattering 
matrix can be generated from Z(J) as follows: 



S p( E p)ab - - t 



dZ(J) 



dJ p (B) t 



L 



dZ{J) 



,7=0 



dJ P (F) ba 



L. 



(3.2) 



,7=0 



where S p (E p ) — diag [S(Ei), £?t(£?a)] an( l L = diag (1,-1). Using standard procedure g], we can represent the 
average of Z( J) as an integral over a 4 x 4 graded matrix field a: 



Z(J)= / d[a] exp{£(<7)}. 



(3.3a) 



where 



C(a) 



N 
2A 5 



trg (cr 2 ) — trgln 



(3.3b) 



Here, trg denotes the graded trace, V(a) = E — a — Ho, co = uo — i0 + , J p (g) = diag [Ji(B), Ji(F), J 2 (B), J 2 (F)], and 
L g = diag (1, -1, 1, -1). We have defined E = (E 1 + E 2 ) /2 and u = E 2 - E x . 
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In the limit N — > oo, this integral can be done with the use of the saddle-point approximation. We are interested 
in correlations involving energy differences ui of the order of the mean level spacing ~ 0(N^ 1 ). Hence, we expand 
C(a) in powers of w: 



N 



C(<r) ~~^ tT S (^ 2 ) -trgln [D(a)] - trgln 



trgln 1 1 - 2tt [V(a) + iirWW^ L] _1 L g WJ p (g)W^ + ^- trg V(a)~^ L 



(3.4) 



I t sh ould be noted such an expansion is not possible for WW^ because W^W ~ 0(1). Of the five terms in expression 
( |3.4|) the last three terms are O(l). The first two terms are O(N) and determine the saddle-point a sp . To derive the 
saddle-point equation we write Ho and a sp in the forms Hq = U^ 1 diag (ei , . . . , ejv) U and a sp = T^ 1 a b ] ^T : where a 8 ^ 
is diagonal and T has the form 



T 



(l + i 12 i 21 ) 1/2 it 12 

-it 21 (l + t 21 t 12 ) 1/2 



The saddle-point equation reads 



sp 



-V 

N ^ 



E-(7r 



(3.5) 



(3.6) 



For ordinary variables (rather than matrices), Eq. ( p.6[ ) has the N + 1 solutions. The N — 1 of which are real, and the 
remaining two may have non-zero imaginary parts according to the values of E. Taking the two complex solut ion s 
(r ± iA) we obtain a s ^ = r — iAL. The explicit expressions of r and A are not available be caus e Eq. ( |3.6| ) 
becomes in general an (JV + l)-th polynomial. Several references discussed the properties of Eq. (3J3) (see, e.g., 
j|,[ll]). Here after we consider the case where A ~ O(l). From the relation between A and the average leve l density 
p (Eq. ( 3.8a )), this means E lies far away from the edge of the spectrum. Substituting a sp for a in Eq. (3^), we find 



C{a sp ) « -trgln [l + inV(a sp y 1 WW'< L 



trg 



V{a sp y 1 L 



-trgln jl - [V{a sp ) +iirWW f L] 1 WL g J p {g)\ . 



The one-point functions p(E) and S(E) are evaluated at the saddle-point. We thus have 



P(E) 



NA 
^A 1 



and 



S P (E) = 1 - 2mW^ [V(a% ) + i*WW*L] WL. 
Using these one-point functions p{E) and S(E) we can write each term of Eq. ( |3.7| ) as follows: 
trgln [l + mV(a sp y 1 WW t L\ = trg In {l - [S P (E) - l] LAI) , 



(3.7) 



(3.8a) 



(3.8b) 



(3.9a) 



CO 



and 



Here. 



trg 



V(a sp ) 1 L\ = -2inLU-p(E) trg (t 12 t 21 ) 



[V(a sp ) + in WW 11 L] 1 W = i {I 1 " 1 !, \S p (E) - l] 1 T-T~ 1 Mt\ 



M 



U 2 t 



12121 



-% t 21 (1 + ti 2 t 21 ) 



1/2 



-i t 12 (1 + £21*12 
—t 2 iti 2 



,1/2 



(3.9b) 



(3.9c) 



(3.10) 
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and we used the property TLT 1 = L + 2M. More explicitly, Eqs. (3.9a) and (3.9c) can be expressed with the use 
of ti2 and t2i as follows: 



R.H.S. of Eq. (3.9a) = trgln (1 + T 12 t 12 t 2 i) 
= trgln (1+2^12*21) 



(3.11a) 



and 



R.H.S. of Eq. (3.9c) = i 



S(E)(l + T 21 t 12 t 21 y 1 - I 

i tai (i + t 12 t 21 ) 1/2 (t 21 - 1 + t 12 t 21 y 



-i t 12 (i + t 2l t 12 y /z (T12 + t 21 t 

1- S^{E)(1 + T 12 t 21 t 12 )- 1 



(3.11b) 



where T12 = 1 - S(E)S^{E) and T 21 = 1 - S*(E)S(E) (T|. (The derivation is given in Appendix.) Thus we find that 
all t he de pendence of Eq. ( |3~7| ) on W and Hq is completely absorbed in p(E) and S(E). Equations (3.9a), ( 3.9tp , 
and (3.9c) show universality. The explicit forms of correlation functions fp are identical with those for the Gaussian 
ensemble and are found in appropriate references (see, e.g., |p1,|l3|). 



IV. SUMMARY 



For the sake of simplicity, we presented the derivation for only the unitary ensemble. In either of the orthogonal and 
the symplectic ensemble, the internal structures of ii2 and £21 differ from the unitary case. However, our derivaton 
is completely independent of such structures and applies equally to the orthogonal and the symplectic ensembles. 
Taking the generating function: 



Z(J) = 



max{m,n) 

n 

9=1 



det [D p (E p )+WJ«(F)Wi] 
det [Dp(E p ) - W4{B)W^] 



(4.1) 



we can show Eq. (1.2) for m > 2 or n > 2 along exactly parallel lines. 

In summary, we have shown that the local universality in the bulk scaling limit still holds for the ^-matrix correlation 
functions even though unitary invariance is broken by the addition of a deterministic matrix to the ensemble. The 
starting random matrix model contains parameters W and Hq which are specific to individual systems. After ensemble 
averaging, these original parameters are compl etely absorbed into much fewer parameters S(E) and p{E). Thus the 
S'-matrix correlaton functions of the type Eq. (1.2) have universal forms which are independent of Hq but for S(E) 
and p(E) and are determined only by the symmetry of the ensemble. This holds for all the three symmetry classes 
(orthogonal, unitary, and symplectic). The derivation can be similarly applied to the spectral correlation functions. 
Thus we have extended the previous results by Brezin and Hikami [|| to the orthogonal and the symplectic ensembles 
though only two-point functions are considered. 

The present results were derived under the restrictions that the correlation functions contain only two values of 
energy, E\ and E 2 , and that V(tp) has a Gaussian form. It is a natural conjecture that the universality of the ^-matrix 
correlation functions holds even if these two restrictions are removed. The increase of the number of energy arguments 
makes the structure of <j s ^ more complicated. With this point taken properly into account, the similar derivation 
is probable. The extension to the general form of V(cp) seems less trivial because we are no longer able to use a 
Hubbard- Stratonovich transformation in order to introduce a graded matrix a. The simlar procedure used in Ref. Q 
may be incorporated into present derivation. 
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APPENDIX: DERIVATION OF EQS. (3.11a) AND (3.11b) 



To derive Eqs. fl3.11a| ) and ( |3.11b| ), we use the fact that for any analytic function F, we have ti 2 F(t 2 iti 2 ) 
F(ti 2 t 2 i)ti 2 . Using the identity 
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a b 



trgln ( ^ , j = trgln (a — bd 1 c) + trgln (d) . 



(Al) 



we obtain Eq. (|3.11a| ). For the abbreviation Y = {t^L [S P (E) - l] 1 T - T^MtX \ using the property T^ 1 
LTL, we get the following equation 



[S(E)-1] l +t 12 t2iA it 12 (l + t 21 ti 2 ) 1/2 A 



it 21 {l + t 12 t 21 ) 1/2 A - S^(E)-1 



where A = {S(E) - l] 1 + &(E) - 1 



a b 
c d 



+ 1. With the use of the formula 

(a-bd- l c)~ l -ar l b(d-ca- l by l 
-d- 1 c(a-bd- 1 C y l (d-ca^by 1 



the [1, 1] block of V can be written as follows 



Y u = A- 1 
Using the property 



5t(£)-l + t 12 t 21 A U [S(E) - 1] '^l- 1 



^12^2 



.4 



St(E) - 1 S't(£;)S , (£;) - 1 [5(S) - 1] 



= [5 , ( J B)-l] S(E)Si(E) - 1 St(£?)_i 
we obtain the [1, 1] block of Eq. ( |3.11b|) . Similarly the [2, 1] block of Y can be written as follows 

,1/2 



Y 2 i = i t 2 i (1 + 1 12*21 J 



St(B)-l 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



Substituting Eq. ( |AJ ) for Fn in Eq. ( |A6|) , we obtain the [2, 1] block of Eq. ( |3.11b[ ). The other blocks are obtained 
along exactly parallel lines. 
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